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TRIPPING  OF  STIFFENED  PLATES  USING  A  REFINED  BEAM  THEORY 


D.A.  Danielson 

Mathematics  Department 

Naval  Postgraduate  School 

Monterey,  CA   9  3  943 


ABSTRACT 

The  subject  of  this  paper  is  the  buckling  behavior  of  thin 
bar  stiff eners  attached  to  plates  which  are  subjected  to  normal 
pressure  loadings.  The  stiff eners  are  mathematically  modeled  by 
a  nonlinear  beam  theory  recently  derived.  An  analytical  solution 
is  obtained  to  the  beam  buckling  equations.  The  curvature  of  the 
centerline  along  the  beam's  base  at  the  buckling  load  is 
expressed  in  terms  of  the  beam's  thickness,  height,  and 
rotational  stiffness.  Analytical  results  are  compared  with  an 
experiment  recently  performed. 


I.    INTRODUCTION:   STATE  OF  THE  ART 

Stiffened  plates  and  shells  are  a  basic  structural  component 
of  submarines,  ships,  and  aircraft.  These  structures  are 
designed  with  generous  safety  margins  against  overall  collapse 
triggered  by  frame  yielding  or  tripping.  Tripping  is  a  lateral 
torsional  buckling  occurring  in  flexurally  stiff  frames  which 
have  low  lateral  rigidity.  The  object  of  analytical  work  is  to 
determine  design  criteria  to  inhibit  tripping  at  any  stress  less 
than  yield.  Tripping  reduces  structural  integrity  and  may 
initiate  failure  of  the  entire  structure  by  general  buckling. 

Surprisingly  little  material  exists  in  the  literature  on  the 
subject  of  the  lateral  instability  of  stiffeners  welded  to 
continuous  plating.  There  are  few  studies  based  on  theories" 
simple  enough  to  have  analytical  solutions.  Earlier  work  is 
summarized  by  Bleich  (1952).  Kennard  (1959)  studied  initially 
curved  stiffeners.  Adamchak  (1979,1982)  pointed  out  the 
importance  of  rotational  constraint  on  the  buckling  load.  Van 
der  Neut  (1982)  developed  a  theory  for  a  Z-stiffened  panel  in 
compression  that  could  be  solved  with  a  pocket  calculator.  More 
accurate  codes  requiring  powerful  computers  were  developed  by 
Smith  (1968)  and  Wittrick  (1968)  based  on  folded  plate  analysis. 
Bushnell  (1985)  also  modeled  the  rings  on  cylindrical  shells  as 
plates. 

In  the  past  there  have  been  few  experiments  in  which 
stiffeners  attached  to  shells  have  been  allowed  to  buckle.  Smith 
(1975)   tested   the   compressive   strength   of   ship   grillages. 


Recently,  at  the  Naval  Postgraduate  School,  a  series  of  stiffened 
plates  have  been  subjected  to  static  or  dynamic  pressure  loads 
sufficient  to  cause  tripping.  Each  plate  was  rectangular  in 
cross  section  and  fixed  at  its  boundaries.  A  single  narrow- 
flanged  stiffener  was  attached  at  its  base  to  the  center  of  each 
plate  and  free  at  its  ends.  Measurements  were  made  of  strains 
and  deflections  versus  pressure,  as  reported  by  Budweg  and  Shin 
(1987)  .  At  low  pressures  a  plate-stif fener  simply  bowed  out 
symmetrically,  but  above  a  critical  buckling  pressure  the 
stiffener  rotated  about  its  base  and  deformed  unsymmetrically 
(see  Fig.  1) . 

In  this  paper  we  mathematically  simulate  these  experiments. 
Our  analysis  is  based  upon  the  following  assumptions: 

(i)     The   stiffener   is  rectangular  in  cross  section  with 

t2     h2 
thickness  t,  height  h,  and  length  i.       -y  and  -^-are 

h 

negligible  compared  to  1. 

(ii)    The  reference  line  along  the  center  of  the  stiffener' s 

base  undergoes  a  vertical  displacement  w  and  negligible 

horizontal  displacement.  YL    is  negligible  compared  to  1. 

h 

(iii)   Each  stiffener  cross  section  does  not  distort  in  its 
plane  and  remains  normal  to  the  reference  line. 

(iv)    The  deformation  normal  to  the  plane  of  a  stiffener  cross 
section  is  egual  to  the  product  of  the  warping  function 
of  the  Saint-Venant  torsion  theory  times  the  twist  of 
the  stiffener. 

(v)     The  stiffener  material  is  linear,  isotropic  and  elastic. 


(vi)    The  plate  does  not  participate  in  the  buckling  of  the 

stif fener. 
(vii)   At  the  buckling  load  the  curvature  of  the  reference  line 
is  a  constant. 
Assumptions   (i)-(v)   are   specializations   of   the   ones 
previously  used  in  deriving  a  refined  nonlinear  beam  theory  used 
to  model  helicopter  rotor  blades  (Danielson  and  Hodges,  1987, 
1988)  .    Hence  we  can  use  that  theory  for  the  present  problem. 
Assumptions  (vi)-(vii)  uncouple  the  beam  problem  from  the  plate 
problem.   Thus  the  mathematical  model  reduces  to  one  dimension, 
in  contrast  to  previous  analyses  in  which  the  stif fener-plate  was 
modeled  by  two-dimensional  plate  theories.   As  a  consequence  of 
the  simpler  formulation,  we  will  be  able  to  obtain  an  analytical- 
solution  to  the  equations. 


II.   NONLINEAR  BEAM  EQUATIONS 

In  this  section  we  reproduce  relevant  formulas  from  our 
previous  papers  (Danielson  &  Hodges,  1987,1988),  applied  to  the 
present  problem.  We  retain  the  same  notation  as  in  these  earlier 
papers. 

The  centerline  along  the  base  of  the  undeformed  beam  is 
called  the  reference  line  r  (refer  to  Fig.  2).  The  Cartesian 
coordinates  of  a  point  in  the  undeformed  beam  are  denoted  by 
(Xi,X2,X3),  where  x^  denotes  distance  measured  along  r  from  the 
middle  of  the  base,  x2  denotes  distance  measured  normal  to  r 
parallel  to  the  plate,  and  x3  denotes  distance  measured  normal  to 
the  plate  from  the  base.   At  each  point  an  orthogonal  reference 

triad  (b^b^b^)  tangent  to  the  coordinate  curves  is  defined," 

r  ... 

with  h±    parallel  to  the  X;l  axis.   The  position  vector  to  points 

in  the  undeformed  beam  is  then  given  by 

r  =  Xib[  (1) 

(The  repeated  index  i  is  summed  from  1  to  3.) 

After  deformation  the  locus  of  material  points  on  the 
reference  axis  is  denoted  by  R.  Now  at  each  point  along  R  define 
an  orthogonal  reference  triad  (b5",b§,b§)  tangent  to  the  deformed 
coordinate  curves.  Also,  let  b  denote  an  intermediate  unit 
vector  in  the  direction  of  the  principal  normal  to  R.  It  follows 
from  assumptions  (i)-(iii)  that 


*-  R   *_ r    .  *_ r 
bi  =  bi  +  w'bo 


r       r 
b  =  b3  -  w'b^ 

R  r  (2) 

b2  =  cos  6   b2  +  sin  9  b 

R  r 

b-i  =  -  sine   bo  +  cos  9  b 


R 
Here  9   is  the  angle  of  rotation  between  b  and  b3 ,  and  primes 

denote  differentiation  with  respect  to  x±.       Expanding  cos  9  and 

sin  9  in  powers  of  9  and  retaining  up  to  quadratic  terms  (higher 

order  terms  are  not  needed  for  the  buckling  analysis) ,  we  obtain 

from  (2) 


r        2 
b2R  =  -w'Bb!  +  (l-^-Jb!  +  9b3 

(3) 
2  2 

bf   =  -(l-^L.)w'bJ  -  0  b£  +  (1-V>b3 


The  curvature  vector  of  the  deformed  beam  is  defined  by 

K  =  <±b\    =   ±to\x    (bf)  '  (4) 

The  components  of  the  curvature  vector  are  obtained  from  (2) -(4) 

K!  =  9» 

K2  =  -w"  +  5^1_  (5) 

K3  =  w  e 


It   follows   from  assumptions   (i)-(iv)   that  the  position 
vector  to  points  in  the  deformed  beam  is  given  by 


R  =  X]_bi  +  wb3  +  x2b2  +  x3b3  +  x2x3e»b1  (6) 


Simplified  expressions  for  the  extensional  strain  Y31  and 
transverse  shear  strains  y12  an<3  Y13  are  derived  in  our  earlier 
papers.  Below  we  reproduce  equations  (9) -(10)  from  Danielson  and 
Hodges  (1988),  specialized  to  the  present  problem: 


Yll  =  En  +  E12*3  "  E13*2  +  ^2  +1*3 

Y12  =  E12  "  ^11*3  (7) 

Y13  =  E13  +  2EH*2 


Ell  =  X3K2  "  X2K3  +  X2X3  6" 


-x^k.  +  x3e'  +  x2x3e'i<3 

E12  =  5 


x0k,  +  xo0'  -  x9x  0'k9 
E13  =  -LI i__ 2_3 2  (8) 

-X2K1  +  X29 '   +  X2X30 'K2 


-x  k   -  x  6*  +  x2x36'K3 

*3   =  2 


Substituting  (5)  and  (8)  into  (7),  and  neglecting  small  and 
higher  order  terms,  we  obtain 

2     2 
Yll  =  _X3W"  +  X2X3  9"  +  — —2 "  X2W"6 

x.w-e2    xlx   w" (6')2 

+  — +  — — 

2  2 

-x^w-e1   x0x2e,eM 

T12  =  "V"  +   2  32  O) 

x„xow"0'   x^xirw-e'e" 

Y         ai      2  3  2  3 

Y13  =  x2  9'  +  5 4 


The   other   strain   components  ( "Y22  'Y  23 » Y33)   turn   out   to   be 

negligible  in  the  strain  energy  expression. 

It  follows  from  assumption  (v)  that  the  strain  energy  of  the 
stiffener  is 


A    t  2      2 

2    2   ho    2y     2y 

"s  =    /    t/   o/  !<yIi  +  rrf +  rrl)dx3  ^2  <*i  do) 

"2    '2 


where  E  is  Young's  modulus  and  v  is  Poisson's  ratio.  Substitut- 
ing (9)  into  (10),  performing  the  x2  and  x3  integrations,  and 
neglecting  small  and  higher  order  terms,  we  obtain 

A 

Ws=     /51W  (11) 


2 
Here  Q  denotes  the  terms  in  Ws  which  are  guadratic  in  0  and 
derivatives: 


8 


its 


Q  = 


A 

2 


E,tV 
4L   18 


^W+^.e.,^ 


(12) 


III.  BUCKLING  EQUATIONS 

We  base  our  buckling  analysis  on  the  energy  criterion  of 
elastic  stability.  This  criterion  and  its  application  are 
explained  by  Timoshenko  and  Gere  (1961),  Danielson  (1974)  and 
Wempner  (1981) . 

The  total  potential  energy  of  the  plate-stif fener 
combination  is  the  functional 

P[U, e]  =  Wp[u]  +  Ws[w»,e]  -  V[u]  (13) 

Here  u  denotes  the  displacement  of  points  on  the  top  surface  of 
the  plate;  note  that  the  vertical  displacement  of  the  plate  along 
the  reference  line  is  equal  to  w.  Wp  denotes  the  strain  energy 
of  the  plate;  Wp  is  a  functional  of  u.  Ws  is  the  strain  energy 
of  the  beam,  given  by  (11) -(12)  ;  Ws  is  a  functional  of  the  beam 
curvature  w"  and  rotation  9  .  V  is  the  work  of  the  external 
loads,  which  is  the  hydrostatic  pressure  times  the  volume  between 
the  undeformed  and  deformed  plate;  V  is  a  functional  of  the  plate 
displacement  field  u  only. 

The  prebucklinq  equilibrium  state  I  in  the  plate-stif fener 
is  denoted  by  (u,w",0=O).  The  potential  energy  in  the 
prebuckling  state  I  is  thus 

PI  =  wptu]  +  ws[w"'°]  "  vtul  (14) 


10 


Invoking  assumption   (vi) ,   we  consider  a  small  deviation 
6^0  from  the  prebuckling  state  of  the  plate-stif fener .    The 
potential  energy  in  this  alternate  state  II  is 

pII  =  wptu]  +  Ws[w»,6]  -  V[u]  (15) 

From  (11),  (1^-)  and  (15),  the  change  in  potential  energy  is  thus 

PlI  "  Pi  =  Q[w",6]  (16) 

According  to  the  energy  criterion  of  elastic  stability,  when  the 
curvature  reaches  a  critical  value  w" ,  there  exists  a  bifurcation 
buckling  mode  0  satisfying 

Q[w",6]  =0,     Q[w",G^]  >  0  (17) 

As  a  consequence  of  (17) ,  the  buckling  mode  is  determined  by  the 
variational  equation 

5Q  =  0  (18) 

Taking  the  variation  of  (12),  and  integrating  by  parts,  we 
obtain 


6Q=  ^e'«e-]  J-  /|[#T,,^^e'  -3T^,]6e,dxi  =  0   (19) 

~2       ~2 


11 


In  order  for  (19)  to  vanish  for  arbitrary  6  0'  ,  0  must  satisfy  the 
differential  equation 


■■■  +  (2hw" L__)e-  =  o  (20) 


2         2 

t    (l+v)li 


and  the  boundary  conditions 


""(f)  =  0"(-  f)  =  0  (21) 


In  order  to  obtain  an  analytical  solution  to  the 
differential  equation  (20) ,  we  invoke  assumption  (vii) .  The 
solution  to  the  eigenvalue  problem  (20) -(21)  is  then 

2nTTx-, 
A  cos[ = ] 

( 

'  or  (22) 


I        (2n+l)7Tx1 
B  sin[ » 1 


where  A  and  B  are  arbitrary  constants  and  n  is  any  positive 
integer.   Furthermore,  the  curvature  is 

„  2  2  2^.2 

—       2t  n  tt  t 

w"  =  j  +  2  (23) 

3(l+v)h  9h£ 


For  small  values  of  n,  the  underlined  term  in  (23)  is 
negligible.  Hence  the  wavelength  of  the  buckling  mode  is  not 
uniquely  determined.  The  critical  curvature  of  the  reference 
line  at  the  buckling  load  is  thus 


12 


2 

3  (1  +  v) 


2t'_ 

3 


w"  =  — T  (24) 


The  maximum  value  of  the  compressive  strain  at  the  buckling 
load  is  obtained  by  setting  9=0  and  x3  =  h  in  the  first  of 
equations  (9) : 

|Yii I      =  h™"  (25) 

1  ' 11 1  max 


13 


IV.   ROTATIONAL  CONSTRAINT 

In  the  preceding  analysis  the  beam  cross  section  was  allowed 
to  rotate  freely  about  the  reference  axis  (simple  support) .  In 
this  section  we  determine  the  effect  of  rotational  constraint  on 
the  buckling  load.  The  method  used  is  explained  by  Timoshenko 
and  Woinowsky-Krieger  (1959)  and  Ugural  (1981) . 

Resistance  to  rotation  can  be  modeled  by  considering  the 
base  of  the  beam  to  be  supported  by  a  foundation,  itself  assumed 
to  experience  elastic  deformation.  We  see  from  (3)  and  (6)  that 
a  point  on  the  base  of  the  beam  undergoes  a  vertical  deflection 
(w  +  X26).  The  foundation  reaction  forces  are  assumed  to  be 
K(w  +  x26).  Here  K  is  a  constant  called  the  modulus  of  the 
foundation  and  has  the  dimensions  of  force  per  unit  area  of  the" 
base  per  unit  deflection.  The  strain  energy  Wf  due  to 
deformation  of  the  foundation  is  then 

1  t 

2  2 

Wf[w,  6]  =  |     /      /   (w  +  x  6)2  dx   dx, 

it  l 


1  1 

■  ¥    /  »2  d*i  +  t£    /  9*  d*i  (26) 

~    2  2 


It  follows  from  the  arguments  in  the  preceding  sections  that 
the  underlined  term  in  (26)  must  be  added  to  the  expression  (12) 
for  the  quadratic  terms  in  the  change  of  potential  energy  Q. 
Then  solving  (17)  for  the  critical  curvature,  we  obtain 


14 


I 

r2   „   ^3,  3 


/{f^e")2  +  3(I^)(e,)2]+^e2^l 


w"  <  Jl _ (27) 

^IT      I2   (e')2  ^1 
~2 
We  assume  that  the  buckling  mode  can  be  approximated  by  our 

previous  solution  (22),  which  is  exact  when  K  =  0: 


2nTTx, 
=  A  sin[ jr-±]  (28) 


where  A  is  an  arbitrary  constant  and  n  is  any  positive  integer. 

Substitution  of  (28)  into  (27)  then  yields  bounds  on  the  critical 
curvature: 

2t2       .„       2t2    ,  4tV   2  +  C*  (29) 


<  r  "  <  = =r  + 


T  ^  w   ^  J  T     2  v       2 

3(l+v)hJ        3(l+v)h     9h£        n 

where 


1    /3K£ 
C  =  — T  W j  (30) 

4tt    V  Eh 


The  best  estimate  for  wH  is  obtained  by  choosing  the  value 
of  n  which  minimizes  the  underlined  term  in  (29).  When  C  <  2, 
the  minimum  value  occurs  for  n  =  1,  and  the  underlined  term  is 
negligible.  When  C  >  2,  the  minimum  value  occurs  for  n  >  1,  but 
is  still  negligible  until  n2  is  large.  When  n2  is  large,  we  can 
treat  n2  as  a  continuous  variable  and  set  the  derivative  with 
respect  to  n2  of  the  underlined  term  egual  to  zero,  which  results 


15 


in  a  value  of  n2  =  C.   Substitution  of  n2  =  C  and  (30)  into  (29) 
then  yields  the  best  bounds 


2t' 


3 (l+v)h" 


<  w   < 


2t' 


3(l+v)h 


2t' 
9h' 


3Kh 


(31) 
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V.    EXPERIMENT 

In  this  section  we  compare  our  formulas  with  an  experiment 
performed  at  the  Naval  Postgraduate  School.  The  experimental 
results  are  taken  from  Budweg  (1986)  and  Budweg  and  Shin  (1987). 

A  rectangular  plate  with  a  narrow-flanged  stiffener  was 
machined  out  of  a  single  large  blank  of  6061-T6  aluminum.  The 
dimensions  of  the  resulting  plate  and  stiffener  are  shown  in  Fig. 
3.  A  strongback  was  bolted  to  the  test  panel.  The  test  panel 
cavity  was  gradually  filled  with  water.  Measurements  of  the 
strains  and  deflections  at  the  bottom  of  the  plate,  and  the 
strain  at  the  top  of  the  stiffener,  were  taken  at  various 
hydrostatic  pressures. 

The  experimenters  judged  that  tripping  of  the  stiffener" 
initiated  at  a  deflection  of  about  three  plate  thicknesses,  when 
the  curvature  of  the  plate  at  the  stiffener  location  and  the 
compressive  strain  at  the  top  of  the  beam  was  approximately  .013. 
The  plate  was  loaded  to  a  maximum  vertical  deflection  of  about 
four  plate  thicknesses.  It  was  observed  that  the  vertical 
deflection  of  the  plate  was  always  symmetric  about  the  plate's 
center  lines.  After  release  of  the  pressure  there  was  permanent 
plastic  deformation  remaining  in  the  plate,  but  no  deformation  of 
the  stiffener  out  of  the  vertical  plane. 

Let  us  now  calculate  the  critical  curvature  obtained  from 
our  analysis  which  ignores  rotational  constraint.  Substituting 
v  =  .33  and  the  dimensions  of  the  stiffener  cross  section  shown 
in  Fig.  3  into  the  formulas  (24)-(25),  we  predict 
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2 

=   .0055      (units  of  inch"1) 


3(1.33)  (1.125)3 


(32) 


Y,  -,  I  =     (1.125)  (.0055)    =    .0062 

1 11 'max 


The  predicted  critical  curvature  and  strain  are  less  than  half  of 

the  measured  values  .013. 

Let   us   examine   the   factors   we   have   neglected   in   our 

analysis,  to  see  which  could  create  a  significant  increase  in  the 

predicted  buckling  load: 

(i)  The  beam  cross  section  was  not  rectangular  and  the 
neutral  axis  in  bending  was  not  at  the  top  of  the  plate. 
However,  fitting  the  cross  section  with  a  more  accurate 
T-shape,  and  assuming  the  beam  is  15/16"  high 
(neglecting  the  fillet  at  the  base  of  the  beam)  ,  we 
obtain  a  critical  curvature  of  only  .0061. 
(ii)  The  prebuckling  state  of  the  stiff ener  was  nonlinear  so 
larger  vertical  and  nonnegligible  horizontal 
displacements  occur.  But  replacing  assumption  (ii)  by 
the  less  stringent  assumption  that  w/h  <  1,  and 
repeating  the  derivation  in  section  III  retaining 
nonlinear  terms  in  w,  we  find  no  substantial  change  in 
(32)  . 
(iii)  The  critical  curvature  w  "  was  not  constant.  However, 
the  experimental  measurements  showed  significant 
variation  in  w"  only  near  the  ends  of  the  stiffener. 
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(iv)  The  structure  had  geometrical  and  material 
imperfections.  But  imperfections  usually  lower  the 
buckling  load. 

(v)  The  beam  cross  section  deforms  in  its  plane.  But 
inplane  extension  and  contraction  have  negligible 
effects  on  the  critical  curvature.  And  allowing  the 
cross  section  to  bend  cannot  increase  the  buckling  load. 

(vi)  The  stiffener  was  partially  restrained  by  the  plate 
against  rotation  at  its  base.  Remembering  the  resux- 
(31)  of  section  IV,  we  must  conclude  that  the  rotational 
restraint  is  large  enough  to  be  a  significant  factor  in 
the  tripping  of  the  tested  plate-stiff ener. 
We  can  estimate  the  magnitude  of  the  rotational  restraint  by 
solving  (31)  for  K: 

K  >   E,9h3w"  -  -I-)2  (33) 

K  -  3h(   2t2    1+V  (33) 

Substituting  w"  =  .013  and  the  beam  dimensions  into  (33),  we 
obtain 

K  >  3E  (34) 
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VI.   CONCLUSION 

As  a  consequence  of  the  theoretical  formula  (31),  we  can 
draw  the  following  conclusions  about  the  tripping  behavior  of 
stiffened  plates: 

1.  The  critical  curvature  "w11  does  not  depend  on  the  length 
i    of  the  beam. 

2.  Beams  with  smaller  thickness  to  height  ratios  t/h  trip  at 
smaller  curvatures. 

3.  The  tripping  point  depends  very  much  on  the  restraining 
stiffness  K. 

Our  mathematical  analysis  has  treated  only  the  stiffener  and 
has  not  considered  its  interaction  with  the  supporting  structure. 
The  analysis  could  be  improved  by  solving  the  prebuckling  problem 
for  the  plate-stif fener,  a  task  which  is  best  done  on  a  computer. 
Additional  kinematical  and  material  nonlinearities  could  be 
included.  The  pressure  which  causes  tripping  could  be  calculated 
and  compared  with  experiment. 

Further  experiments  need  to  be  done.  It  was  difficult  to 
determine  when  and  if  tripping  occurred  from  the  measurements 
that  were  made.  Future  experiments  should  allow  visual 
inspection  of  a  stiffener  while  it  is  buckled.  Additional 
experimental  measurements  could  be  made  to  determine  the 
rotational  stiffness. 

The  advantage  of  the  present  methods  is  that  they  lead  to 
simple  analytical  formulas  which  reveal  the  dependence  of  the 
buckling   point   on   geometrical   parameters.     It   would   be 
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interesting  to  see  if  these  analytical  techniques  could  be 
applied  to  other  problems,  such  as  the  buckling  of  initially 
curved  shells  with  stiffeners  of  T  or  Z  cross  sections. 
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FIGURE  CAPTIONS 

Figure  1:   Stiff ener  shapes  when  plate  is  under  various  pressures 
Figure  2:   Geometry  used  in  mathematical  analysis 
Figure  3:   Dimensions  of  experimental  model 
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